Asymmetric radiating brane-world 
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At high energies on a cosmological brane of Randall-Sundrum type, particle interactions can 
produce gravitons that are emitted into the bulk and that can feed a bulk black hole. We generalize 
previous investigations of such radiating brane-worlds by allowing for a breaking of ^-symmetry, 
via different bulk cosmological constants and different initial black hole masses on either side of the 
brane. One of the notable features of asymmetry is a suppression of the asymptotic level of dark 
radiation, which means that nucleosynthesis constraints are easier to satisfy. There are also models 
where the radiation escapes to infinity on one or both sides, rather than falling into a black hole, 
but these models can have negative energy density on the brane. 



I. INTRODUCTION 

Brane-world models of the universe allow us to ex- 
plore the cosmological implications of extra-dimensional 
gravity. In these models, matter is confined to a four- 
dimensional "brane" while the gravitational field propa- 
gates in a higher-dimensional "bulk" spacetime. Of par- 
ticular interest are the Randall-Sundrum (RS) type cos- 
mological models with one extra spatial dimension 
that is infinite. At low energies the 5D graviton is local- 
ized at the brane, because of the strong curvature of the 
bulk due to a negative bulk cosmological constant, and 
general relativity is recovered with small corrections. At 
high energies, gravity behaves as a 5D field and there 
are large deviations from general relativity. The effec- 
tive 4D Planck scale on the brane, M4 « 10 19 GeV, 
can be much higher than the true 5D Planck scale M5 
(10 8 GeV < M 5 < M^. 

The background RS-type cosmology is a 
Schwarzschild-anti de Sitter (AdS) bulk with a 
Friedman-Robertson- Walker brane that is moving 
away from the black hole as the universe expands. The 
mass parameter m of the bulk black hole produces a 
Coulomb effect at the brane. This Coulomb effect is 
felt on the brane as a "dark radiation" term m/a 4 in 
the modified Friedman equation on the brane. One 
of the corrections to the background model that can 
be considered, is to incorporate the production of 5D 
gravitons in particle interactions 0. This can become 
significant at very high energies. The gravitons are 
emitted into the bulk and feed the black hole, so that 
771 > 0. At low energies the effect is negligible and m 
asymptotes to a constant value. The general problem 
is highly complicated, but a simplified model can be 
constructed by neglecting the back-reaction of the 
gravitons on the bulk geometry, and by assuming that 
all gravitons are emitted radially 3] (some progress in 
the case with non-radial emission has been made Q). 
This leads to a Vaidya-AdS bulk (previously considered 
as a model for thermodynamic interactions between the 
brane and the black hole 5j). The dynamical equations 
for the simple radiating brane-world may be solved 
exactly The black hole mass grows monotonically 



and tends to a late-time constant. An expanding brane 
is always outside the black hole horizon Q. 

Here we consider a further generalization of the radi- 
ating brane-world model, by allowing for asymmetry in 
the bulk, i.e., different cosmological constants and initial 
black hole masses on either side of the brane. Such asym- 
metry has been previo usly investigated for non-radiating 
brane-worlds [1HHII1H0- ( The equations for the 
completely general case are derived in [l4|.) We find that 
for radiating brane-worlds, the introduction of asymme- 
try leads to some interesting new features, in addition 
to those already identified previously for non-radiating 
models. 



II. THE BULK AND BRANE 

The bulk metric is Vaidya-anti de Sitter [jj, 

ds 2 = -f(v, r)dv 2 + 2e dvdr + r 2 dx 2 , (I) 

where v = const are plane-fronted null rays which are 
ingoing for e = 1 and outgoing for e = — 1, and 



f(v,r) 
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(2) 



Here A is the bulk cosmological constant, and m is the 
mass parameter of the black hole at r = 0. We assume 
that the brane is outside the horizon (/ > 0), which 
requires negative A for non-negative m. The metric is a 
solution of the bulk field equations 



^G ab + A(% ab = K 2 &T a 



(3) 



where k 2 = 8tt /-Mf, and the bulk stress tensor has null- 
radiation form, 



( 5 >T ab = iP(v,r)k a k b . 



(4) 



Here ip gives the power density in the null radiation that 
flows along the null rays generated by k a , where k a oc 6 a v . 
The field equations require that 



3t3 2 dm 
r 3 dv 



(5) 



2 



where v = dv/dt and t is cosmic proper time on the 
brane. The brane trajectory is r — a(t), and we have 
taken the brane to be spatially flat. The ray vector has 
been normalized so that k a u a = —1, where u a is the 
brane 4- velocity. From u a u a = — 1 we obtain 



fv = er + 



1/2 



(6) 



Up to this point everything applies for either side of the 
brane, i.e. for A = A±,m — m±,e = e±. 

At high energies on the brane, p > A, where A is 
the brane tension, particle interactions produce gravi- 
tons that escape into the bulk. Assuming that all such 
gravitons are emitted radially, and that we can neglect 
their back-reaction on the bulk geometry, we can use the 
Vaidya-AdS metric to model the bulk. (Corrections in- 
troduced by allowing for non-radial emission require nu- 
merical computation 4].) The function ijj is determined 
in the radiation era (p = p/3) by kinetic theory argu- 
ments as 0: 



K 2 

^ap 2 . 
12 P 



(7) 



where a is a small dimensionless constant. Because of 
graviton emission, the brane energy is not conserved, 



p + AHp = -2i/> , 



(8) 



where H — a /a is the Hubble rate. 

The radiation leaving the brane is the same on both 
sides of the brane, but the different curvature on either 
side affects the propagation of the radiation. The pa- 
rameters e± describe where the radiation is going: e = 1 
indicates radiation into the black hole, e = — 1 indicates 
radiation to infinity. There are four cases: 

• For e± = 1, the brane radiates into black holes on 
either side of the brane. 

• For e± = —1, the radiation is escaping to infinity 
on both sides of the brane. 

• For the remaining two cases, with e+e_ = — 1, radi- 
ation falls into a black hole on one side and escapes 
to infinity on the other. 

Note that in the non-radiating case (m = 0), there are 
also four cases with the above values of e±, with e = 1 
denoting the half-bulk that includes the black hole (or 
that includes the AdS Cauchy horizons in the case m = 
0), while e = — 1 denotes the half-bulk on the other side 
of the brane. 

These four cases show important differences in the re- 
lation between the brane energy density and the bulk 
geometry. The key equation is Eq. (77) in 0], which in 
our scenario becomes |l5| 



The simplest generalization of the ^-symmetric mod- 
els in 3] is the case with e± = 1, and in this case, Eq. © 
shows that p + A > 0. By contrast, the e± = —1 models, 
in which the radiation escapes to infinity on both sides, 
have p + A < 0. It is not possible to have radiation es- 
caping to infinity on both sides without having negative 
total energy density on the brane. This negativity also 
holds in the non-radiating case and in the ^-symmetric 
case; it is not the asymmetry or the radiation escaping 
that causes the problem, but the fact that the bulk cur- 
vature is decreasing as one moves away from the brane. 
The models with e+e_ = —1 will in general have p + A 
changing sign over time. In summary, 

e± = 1 => p + A>0, (10) 
e± = -1 => p + A < 0, (11) 
e + e_ = — 1 p + A can change sign. (12) 

Equations (j^J , © and JJJ show how the mass function 
evolution at the brane depends on the bulk geometry on 
either side of the brane: 



m± 



"36 
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-H + e ± 



1/2 



(13) 



The average, m = (m + + m_)/2, gives 



30 



-m + a = e + [ a 2 + f + -J- 



1/2 



< (« 2 + /-^) • n-i.) 



1/2 



where 



- 2m A , A 2Ato AA , . , 

/ = 2 7f a > f = 2 7T a ■ 15 

a z 6 a z 6 

The difference, Am = rn + — m_, follows from the square 
of Eq. C3J: 



36 • \ . . niaa 2 p 2 . „ . . 
m + H ) Am = A/ . (10) 
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The generalized Friedmann equation for an asymmetric 
brane is given in |l4j. In our case, 



2_4 ( x) 2 f 9 (A/) 2 

H "~36 {P+ } ~^ + 44a 4 (p + A) 2 



(17) 



This can be employed to eliminate a 2 + f from Eq. 114fl , 
giving 



30 
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m + H = € 



4 (P + A) 



+ Ae 



3A/ 



2k§ (p + A) a 2 



(18) 
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(p + A) = e+^TU + e- V^+f- 



From now on we will assume that p > and A > 0, 
(9) and confine our investigation to the e± = 1 scenario. In 



3 



this case, the projected field equations on the brane p) 
show that the effective gravitational constant is 



K 2 = 



8tt 

m! 



4\ 

6 



(19) 



We also fine-tune the cosmological constant on the brane 
to zero, 



k X + A = . 



(20) 



Then the generalized Friedmann equation (|17|l and 
generalized Raychaudhuri equation jl4j become 



H + H< = -- P 



2m 



9 (A/) 2 



a 4 4 K 4 a 4 (p + A) 2 
(3 + a)p 



2A 



2m 



3AA A/ 



9(p-3A) (AfY 



2 K |a 2 (p + A) 4k|o 4 (p + A)" 



(21) 



(22) 



In the symmetric case we recover the equations of . 
The exact solution is given by [||: 

a = a^t-to^+^/^+^it-h) 1 /^, (23) 

9 = K?A{t-t )(t-h) ' (24) 
a 4 L4(A - 8)(* - ti) + 4^{A - 2(4 + a)}] 

777 = — ■ 

8A 2 (t - i ) 4/(4+a) (* - ti)°/( 4+a ) 
where the constants are related by 



,(25) 



(26) 



In the physically relevant case (i.e., H > 0,m > 0), we 
have t > ti > t and A > 8 + 2a. 



III. SOLUTIONS 

For the general asymmetric case, we define dimension- 
less quantities as follows: 



t 



P = 



-t, H 



H. 



rn 



6 
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or 



(27) 



(28) 



(29) 



The first four variables are the same as those used in 0, 
6], while q and the parameter £ encode the asymmetry. 

We constrain £ by the requirement that A± < 0. This 
means that AA < A, so that the cosmological constants 
on both sides of the brane are negative. Then from 
Eqs. igDJ and 0, we find that 



0<£< 



(30) 



Using these variables, the dynamics on the brane in 
the case e± = 1 are governed by the system of first order 
linear differential equations, Eqs. JEJ, fTH|l . ltT%|) and lf2*2*| . 
In normalized form this system becomes 



P = 



rn 



Q 

H' 



-4Hp-ap 2 , 

1 - H 
p 2 (q + 3^) 



aa p \ p 



(P + 1) 



-H 2 



2p 
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(3 + a) 



2 m 



3£(g + 3£a 4 ) , (p - 3)(g + 3£a 4 ) 
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(31) 
(32) 

(33) 



(34) 



The algebraic constraint from the generalized Friedmann 
equation, 



if 2 



2m 

P(2 + P) + -T 



3^y 



a 4 4a 4 Q5+l) 2 ' 



(35) 



is used to monitor numerical errors. Results from a nu- 
merical integration of this system with values of £ in the 
range of Eq. Ij30|l and a variety of initial conditions for q 
are shown in Figs. 1-5. 




FIG. 1: The scale factor a(t) with q(0) = 0. The curves from 
bottom to top have £ increasing from to | . 

The first three graphs show the effect of increasing £ 
while keeping q(0) fixed. The effect on the scale factor 
is to introduce acceleration at late time. This effect is 
independent of the Vaidya radiation and is also found in 
the non-radiating case Hi. (See 1MII3 

for asymme- 
try in bulk cosmological constants, |l2j for asymmetry in 
(constant) bulk black hole masses, and 0, Qj| for both 
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types of asymmetry). The effects of the radiating brane 
occur predominantly at early times. 

There is also a marked effect on the average mass func- 
tion rh. The increase of H when £ is increased in Eq. Ij35(l 
leads to a decrease in rh' by Eq. (|32|) . which results in 
a suppression of the asymptotic value of rh. Since the 
asymptotic value determines the dark radiation at nu- 
cleosnthesis (where p <C A) .this suppression of rh means 
that nucleosynthesis limits [3j can be more easily satisfied 
when there is asymmetry. 

Figure shows that even if there is no initial differ- 
ence in the black hole mass functions on either side of 
the brane, and given that the radiation from the brane 
is symmetric, a difference in the mass functions will still 
develop during the evolution, because of the different cur- 
vatures on either side of the brane. 

Figures 01 and [SI demonstrate the effect of varying q(0) 
for a fixed £ . Although there is little discernible difference 
in the scale factor for the different values of q(0), the 
increase in H for larger q(0) is still enough to suppress 
the final value of rh. (The shape of the q curve is the 
same as in Fig. |3 for all the initial values.) 

If we assume that the asymmetry is small, i.e. q/a 4 and 
£ are small, then we can perturb away from the symmet- 
ric solution, Eqs. Q23 [l -I|25| ) . and evaluate the early- and 
late-time solutions. At early times, 



70i 
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FIG. 4: The scale factor a(t) with f = \ and q(0) = 0, 10, 
20, 30, 50 and 100, from bottom to top curve. Variation in 
q(0) has little effect on the evolution of a. 



is confirmed in Fig. EI At late times, 



3a 



*2 

a \ t 



(36) 



Be- c / t2 , 



(37) 



where i2 is a constant of integration which must be of 
a suitable size to keep the whole expression within our 
smallness assumptions. This type of decay at early times 



where B, C are constants. The asymptotic value, q 
is evident in Fig. |3 
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FIG. 5: The average black hole mass function m(t) with £ = 
and (curves from top to bottom) 5(0) = 0, 10, 20, 30 , 50 and 
100. 



symmetry. The radiation from the brane propagates in 
different curvature regions on either side, given the dif- 
ference in bulk cosmological constants and initial black 
hole masses on either side. There are four cases, depend- 
ing on whether the radiation falls into a black hole or 
escapes to infinity on either side of the brane. 

We performed numerical integrations of the equations 
in the case where radiation feeds black holes on both 
sides of the brane. One of the key features that emerges 
is the suppression of the late-time dark radiation (the 
asymptotic value of the average mass parameter) due to 
asymmetry. 

These solutions correspond to the case e± = 1 in 
Eq. ©, or e = +1, Ae = in Eq. JTHJ}. The other three 
cases have radiation escaping to infinity on at least one 
side of the brane, but in these cases, the total brane en- 
ergy density p + X is negative or can change sign. Further 
investigation is needed to determine whether any of these 
cases admit physically reasonable models. 



IV. CONCLUSION 

We have derived the dynamical equations for a radi- 
ating brane in a Vaidya-AdS bulk, without assuming Z2 
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